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ABSTRACT 



Aims. We attempt to detect starlight reflected from a hot Jupiter, orbiting the main-sequence star HD 75289Ab. We report a revised 

analysis of observations of this planetary system presented previously by another research group. 

Methods. We analyse high-precision, high-resolution spectra, collected over four nights using UVES at the VLT/UT2, by way of data 

synthesis. We try to interpret our data using diiferent atmospheric models for hot Jupiters. 

Results. We do not find any evidence for reflected light, and, therefore, establish revised upper limits to the planet-to-star flux ratio 

at the 99.9% significance level. At high orbital inclinations, where the best sensitivity is attained, we can limit the relative reflected 

radiation to be less than e = 6.7 x 10"'' assuming a grey albedo, and e = 8.3 x 10"^ assuming an Class IV function, respectively. 

This implies a geometric albedo smaller than p = 0.46 and p = 0.57, for the grey albedo and the Class IV albedo shape, respectively, 

assuming a planetary radius of 1.2 Rjup. 

Key words. Methods: data analysis - Techniques: radial velocities - Stars: individual: HD 75289A - planetary systems 
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1. Introduction 

Since the detection of the first exoplanet orbiting a solar-type 
star, more than 270 exoplanets have been detected. So-called hot 
Jupiters - giant planets only a few solar radii away from their host 
stars - provide the opportunity to detect starlight reflected from 
these planets. Five extended campaigns for the search for re- 
flected light were completed by different groups, all resulting in 
non-detections (e.g. Leigh et al. 2003b, and references therein). 
Upper limits to the planet-to-star flux ratio and to the geomet- 
rical albedo of these planets were established, which provided 
important constraints to models of the planetary atmospheres by 
Sudarsky et al. (2000, 2003). As a result, models that predicted a 
high reflectivity for the planetary atmosphere could be ruled out 
for some of the studied planets. 

Here, we present a reanalysis of observations of the plane- 
tary system of HD 75289A, conducted by Leigh et al. (2003a), 
over four nights in January 2003, using the UV- Visual Echelle 
Spectrograph (UVES) mounted on the VLT/UT2 at CeiTO 
Paranal in Chile. These authors attempted to detect starlight re- 
flected from the hot Jupiter, but were unable to find evidence for 
the planetary signal: they placed a 99.9% confidence upper limit 
on the planet-to-star flux ratio of 4.17 x 10"^, for the spectral 
range A - 402.33 to 522.13 nm, and for an orbital inclination 
i = 60°. 

We noticed, however, that this upper limit was based on er- 
roneous orbital phase information for the planet. In this article, 
we coiTect the upper limit to the planet-to-star flux ratio de- 
termined by Leigh et al. (2003), using using our implementa- 
tion of the modeling approach introduced by Charbonneau et al. 
(1999). Section 2 describes the basic ideas of the search for re- 
flected light. Section 3 provides a brief overview of the science 



data, while Section 4 provides a detailed description of sophis- 
ticated data processing implemented by using the data modeling 
approach. Finally, in Section 5 we present our corrected upper 
limits to the planet-to-star flux ratio. 

2. Reflected light 

2.1. Photometric Variations 

For exoplanets, the enormous brightness contrast between the 
star and the planet constitutes considerable challenge when at- 
tempting to observe some kind of direct signal from the planet. 
For close-in planets such as hot Jupiters, the main contribution 
to the optical flux originates in the reflected starlight and not 
the intrinsic luminosity (Seager et al. 2000). High-resolution 
spectroscopy in the optical utilises the fact that the observed 
spectrum reflected from the planet is essentially a copy of the 
rich stellar absorption-line spectrum. Basically, this spectrum 
is shifted in wavelength according to the orbital radial velocity 
(RV) of the planet and scaled down in brightness by a factor of a 
few times 10"* for hot Jupiters. According to Charbonneau et al. 
(1999), the amount of starlight reflected from a planet which is 
fully illuminated can be described by 



eiA)^p{A)\^ 
a 



(1) 
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where p(A) denotes the albedo of the planet as a function of the 
wavelength A, Rp the planetary radius and a the star-planet sepa- 
ration. Figure[T]shows different albedo spectra, derived using dif- 
ferent planetary-atmosphere models from Sudarsky et al. (2000). 
The planetary radius of HD 75289Ab is unknown; it can, how- 
ever, be estimated from the transiting planets, which provide an 
exact determination of their masses and radii (see Section 2.3). 
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Table 1. Parameters of the star HD 75289A and its planetary 

system 



Fig. 1. Different albedo spectra of atmospheric models (taken 
from Sudarsky et al. 2000) are shown. The irradiated (dots) and 
isolated (dashed) Class IV models describe atmospheres of plan- 
ets with temperatures Teff ~ 1300 K. In contrast to the isolated 
model, the irradiated model assumes that no reflective clouds ex- 
ist in the upper layer of the planetary atmosphere, which results 
in a very low albedo. The solid line depicts a Class V roaster, 
describing the atmosphere of a planet with Teff > 1500 K, hav- 
ing a high reflective silicate cloud deck in the upper layers of the 
atmosphere. 

The orbital radius can be tightly constrained using Kepler's third 
law. 

In most cases, the planet does not appear to be fully illu- 
minated. Consequently, the observed reflected light is reduced, 
depending on the model describing the scattering behaviour of 
the atmosphere, its orbital inclination / e [0°, 90°] and the orbital 
phase (p e [0, 1] of the planet. We note that we adopt the conven- 
tion that = represents inferior conjunction of the planet (for 
i = 90°, it would be the transit position). We apply an empirical 
scattering model of the atmospheres of Jupiter and Venus (Hilton 
1992), which can be approximated by 



where 

l\.2,a\ /l.Sof 

^(q')=0.09 +2.39 

and the phase angle a: 
cos a = - sin / cos 2n^. 
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(2) 
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Finally, the flux of the reflected light from the planet at the orbital 
phase <p follows from equations [1] and |2l 



/(0,/,^) = 6(^)M0,O- 



2.2. Doppler shifts 



(5) 



The planet orbiting its host star produces not only a flux varia- 
tion (Equation|5]l, but also a Doppler shift of the stellar spectrum 
reflected from the planet. The RV semi-amplitude Kp of that shift 
depends on the orbital inclination /, which is unknown for non- 
transiting planets. K^ can be expressed by 



Parameter 


Value 


Error 


Ref. 


Star: 








Spectral type 


GOV 




G89 


V (mag) 


6.35 




VF05 


d(pc) 


28.94 


0.47 


VF05 


M, (Me) 


1.23 


0.10 


VF05 


R. (Ro) 


1.249 


0.022 


VF05 


/'rot (d) 


16.0 


3.0 


UOO 


V sin ;■ (km s ' ) 


4.14 




B06 


Age (Gyr) 


2.2 




VF05 


Planet: 








MpSini(Mj„p) 


0.467 


0.041 


B06 


a(AU) 


0.0482 


0.0028 


B06 


e 


0.034 


0.029 


B06 


K,(kms-') 


0.0549 


0.0019 


B06 


Orbital period (d) 


3.509267 


0.000064 


B06 


r^=o (JD) 


2 450 829.872 


0.038 


B06 



Kp — Ks 



M, 



Mp sin/ 



■ smi. 



(6) 



References: B06 = Butler et al. (2006) and references therein, 
G89 = Gratton et al. (1989), UOO = Udry et al. (2000), VF05 = Valenti 
& Fischer (2005). 



where K^ is the RV semi-amplitude of the star, and Ms and 
Mp sin / are the stellar mass and the minimum mass of the planet, 
respectively. The largest possible amplitude A'p.max occurs at the 
orbital inclination / - 90°. 

The instantaneous RV shift of the planetary signal with re- 
spect to the star depends on the orbital phase <p, 

Vp^KpSm27T<p. (7) 

2.3. HD 75289 A and its planet 

High-precision RV measurements revealed the existence of a 
hot Jupiter orbiting the GO-type main-sequence star HD 75289A 
(Udry et al. 2000). Table [T] summarises the parameters of the 
planet and its host star. We note in passing that the system 
also contains a faint low-mass stellar component, separated by 
*: 621 AU from the primary (Mugrauer et al. 2004). 

The planetary system of HD 75289A seemed to bear the op- 
portunity to detect the reflected starlight for the following rea- 
sons: First, the brightness of HD 75289A in the visual {V = 6.35 
mag) enables a large amount of high-S/N spectra to be acquired 
within a short period of time, which helps to reduce photon 
noise, the dominant noise source. Second, the amount of starlight 
received from this planet is likely to be high due to the extremely 
small distance to its host star of only a = 0.048 AU (cf. Equation 
[U. Third, a GO V star shows a rich absorption-line spectrum 
at optical wavelengths, a prerequisite for the data-synthesis ap- 
proach outlined in Section 4. 

Using data of the stellar mass Ms, the planetary minimum 
mass Mp sin /, and the RV semi-amplitude of the reflex motion 
of the star K^, the maximum possible RV semi-amplitude of 
the planet can be determined (from Equation |6]l to be Kpmn^ = 
148.9 + 15.8 km s"'. An estimate of the planetary radius Rp can 
be obtained from a comparison of the radii determined for other 
hot Jupiters with the transit method. To this end, we selected a 
subsample of known transiting hot Jupiters (Burrows et al. 2007, 
and references therein), based on the following criteria: (i) Planet 
mass between the minimum mass of HD 75289Ab and 2.3 times 
this value (for random orientation of the planetary orbit there is 
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a 90 % confidence that the true mass of HD 75289Ab lies in this 
interval), (ii) Orbital radii similar to that of HD 75289Ab within 
20 %. (iii) Host star spectral types close to GO V (F8 V - Gl V). 
The resulting group of objects contains OGLE 10b, OGLE 1 1 lb, 
HAT-P-lb, XO-lb, HD 209458b. From this sample, we esti- 
mated the planetary radius of HD 75289Ab as the average of 
the radii of these five hot Jupiters finding Rp = 1.2 + 0.14 Rjup- 
The uncertainty corresponds to the scatter in the above sample 
of transit radii. 

3. Data 

We retrieved 684 high-resolution and high-S/N spectra files from 
the ESO science archive, along with the calibration files. The 
observations were conducted over four nights by Leigh et al. 
(2003 a) with UVES mounted on the VLTAJT2. The blue arm 
of the spectrograph was used with the single EEV CCD-44 ar- 
ray detector centered at wavelength A = 475.8 nm providing 
25 full orders and covering the wavelength range A = 402.3 
to 522.1 nm. The integration times were adapted to the seeing 
such that a count rate not exceeding 40 000 ADU per pixel was 
achieved, insuring high S/N levels and at the same time staying 
clear of the saturation level. The effective resolving power was 
R - 43 000. More details of the observations can be found in 
Leigh et al. (2003a). 

These authors selected the observing dates for the best visi- 
bility of the planet, which is in the orbital phase range = 0.30 
to 0.45 and 0.55 to 0.70 (close to superior conjunction), but ex- 
cluding phases about 0.50, where the absorption-line systems of 
the star and the planet blend. However, we found that a one- 
day error in the date used for the phase calculation must have 
been made. This resulted in the fact that the planet was not ob- 
served at the optimum phases. In particular, during two nights 
the observations were taken close to the inferior conjunction 
of the planet, where it is hardly visible. Furthermore, the same 
erroneous phases were used in the subsequent data analysis. 
Consequently, one can expect the corrected upper limit to the 
planet-to-star flux ratio to be smaller than the value determined 
by Leigh et al. Table |2] shows the journal of observations with 
the corresponding phases of the planet (erroneous and corrected 
ones). 

3.1. Data reduction 

The data were reduced using a UVES data-reduction pipeline 
developed by ourselves. For high-quality flat-field correction, 
228 and 369 flat field exposures were combined for the first two 
nights and for the last two ones, respectively. Furthermore, 90 
bias exposures were combined. We retrieved 25 orders of 3000 
pixels from each echelle spectrum. No order merging was ap- 
plied. Furthermore, we identified cosmic -ray hits by way of the 
following procedure: for each spectrum, we compared the flux 
in every pixel with the median flux of the same pixel in the 
three predecessor and the three successor spectra, which had 
been scaled to the same flux as the spectrum under considera- 
tion. We flagged those pixels where the diff'erence exceeded 6cr 
as cosmic-ray hits. These pixels were then excluded from further 
analysis. 

We discarded the most weakly exposed regions of each 
echelle order (the first 300 pixels as well as the last 100 pixels). 
To speed up the data analysis, we co-added spectra into sufli- 
ciently narrow phase bins such that phase smearing of the stellar 
absorption lines originating from stellar RV variation, sub-pixel 
shifts (see Section |4m i, and the barycentric motion of the Earth, 



remained below 0.04 km s ' . An additional criterion for the size 
of the phase bins was that the unseen planetary lines did not 
suffer from smearing in excess of 2 km s"'. In these ways, we 
reduced the number of spectra from 684 to 210. 



4. Data analysis: The data synthesis method 

We model the starlight reflected from the planet as a copy of 
the stellar spectrum, strongly scaled down in brightness and 
Doppler-shifted according to the orbital motion of the planet. 
The HD 75289A spectra have an average S/N of 300 to 600 per 
dispersion element. With the expectation the planet-to-star flux 
ratios of the order of a few times lO"*", it is clear that the re- 
flected spectrum from the planet is deeply buried in the noise of 
the stellar spectrum. The weak planetary signal is increased by 
the large number of spectra, and more importantly, by the com- 
bination of the approximately 1500 absorption lines, achieved 
using the data-synthesis method described below. 



4.1. Step 1: Construction of ttie superspectrum 

To create a high S/N, virtually planet-free superspectrum, we 
first create copies of the original, unmodified object spectra. 
These copies are then co-aligned with respect to the barycen- 
tric velocity of the Earth and the radial velocity of the star, and 
finally summed up. Here, it is important that the observing dates 
cover different orbital phases well distributed such that the plan- 
etary signal in the superspectrum is washed out. In the case of 
the HD 75289A data, nearly half of the observations were con- 
ducted in the orbital phase range = 0.97 to 0.07 (Table |2]i, 
where the planet was close to inferior conjunction and therefore 
at its faintest (e.g. for an orbital inclination of / = 60°, the planet 
appears about 12 times fainter at the inferior conjunction, rel- 
ative to its maximum brightness). Consequently, we used only 
these "planet-free" spectra to create the superspectrum. 

We found that imperfections in the dispersion solution (sub- 
pixel shifts and stretches/contractions) of the different spectra 
led to broadening of absorption lines in the superspectrum. We 
assume that these shifts originate mostly from errors of the wave- 
length calibration and probably to some extent from guiding er- 
rors and variations in the instrumental profile. Consequently, we 
correct each spectrum for these wavelength errors before pro- 
ducing a new version of the superspectrum. To this end, we first 
divide each order of the modified copies of the object spectra into 
chunks of size A^ pixel, and determine the central wavelength Aj 
of each chunk /. For the data analysis of HD 75289A, we set 
A^ - 120, which creates 22 chunks per order By means of the 
Brent algorithm (Press et al. 1992), we determine for each chunk 
the sub-pixel shift si between the superspectrum and the modi- 
fied copy of the object spectrum under consideration. Then, we 
shift each pixel of the modified copy by a value determined by 
calculating a spline through all pairs of Aj, Sj per spectral or- 
der Finally, we add up the re-adjusted co-aligned versions of the 
copies of the object spectra and get an improved version of the 
superspectrum. 

4.2. Step 2a: Modeling tlie stellar signal ... 

Each of the original, unmodified object spectra is modeled using 
the improved superspectrum. In this initial step, we ignore the 
presence of the faint planetary signal in the data. First the super- 
spectrum (model) is corrected for a general linear trend in flux. 
Second, the model is shifted according to the barycentric veloc- 
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Table 2. Journal of observations. The UTC times and the orbital phases of the planets are shown. Note that 0Leigh indicates the 
erroneous phases applied by Leigh et al. (2003), while (^coir refers to the corrected ones. The number of spectra obtained is given in 
the last column. Ephemerides used are; orbital period - 3.509267 d, T^^o (BJD) - 2 450 829.872 (Butler et al. 2006, and references 
therein). 



Night 


UTC start 


(^Lcish 


(peon 


UTC end 


^Lcish 


0corr 


^* spectra 


1 


2003/01/14 01:00 


0.39 


0.68 


2003/01/14 09:30 


0.50 


0.79 


188 


2 


2003/01/15 00:57 


0.68 


0.97 


2003/01/15 09:22 


0.78 


0.07 


173 


3 


2003/01/21 02:40 


0.41 


0.70 


2003/01/21 09:35 


0.49 


0.78 


183 


4 


2003/01/22 03:36 


0.71 


0.00 


2003/01/22 09:41 


0.78 


0.07 


140 



ity of the Earth, the radial velocity of the star, and the aforemen- 
tioned shifts and stretches/contractions in the sub-pixel regime, 
so that the positions of absorption lines of the object spectrum 
and the model are matched. This is achieved using a chunk- 
Brent-spline-approach similar to that described in the previous 
paragraph. We note that in all analysis steps, modifications are 
exclusively applied to the model, but the object spectrum are 
used in their original version, i.e. the data are unchanged. 

At the third stage, the model is scaled chunkwise with re- 
spect to the object spectrum. As we now compare the scaled 
model with each object spectrum, we notice that the widths and 
depths of the absorption lines differ slightly. These differences 
originate most likely in the aforementioned effects of residuals 
wavelength calibration errors, guiding errors, and variations in 
the instrumental profile, and can be corrected by adding a scaled 
version of the second derivative of the object spectrum to the 
model. The scaling factor is determined via x^ minimisation 
(fourth stage). After this, the model is renormalised. The final 
stage is to iterate twice over all these four processes to improve 
the model describing the stellar spectrum. 

4.3. Step 2b: ... and the planetary signal 

For the model of the planetary signal, we use a copy of the im- 
proved model of the stellar spectrum, but scaled down by the 
factors e(A) ^{^, i) and shifted by velocity V^iK^, (p) with re- 
spect to the stellar spectrum. Hence, the two free parameters are 
the planet-to-star flux ratio for the fully-illuminated planet €(A), 
and the orbital inclination /, which corresponds to the RV semi- 
amplitude of the planet K^ - /Tpmax sin / = 148.9 sin / km s"' . 

4.4. Step 3: Evaluation of the models 

We are now ready to add this planetary signal to the improved 
model T of the stellar spectrum and consequently construct the 
model M describing the spectrum of the star plus the reflected 
one from the planet. For each pixel k, M is given by 



Mk^ 



TkJAi) + e(Ak) fi{<P, i) niAk [1 + Vp(Kp, <p) c-']} 
1 + e(Ak) ju(0, /) 



(8) 



where c denotes the speed of light. Varying Kp and e(A), we fi- 
nally search for the best-fit model M to all the object spectra by 
X^ minimisation. The search range for the RV semi-amplitude 
comprised A'p - 40 to 180 km s"' (corresponding to orbital incli- 
nations / - 15° to 90°, plus twice the error of /Tp maxJ see Table 1) 
with a step width of 3 km s"'. This was a good compromise be- 
tween computing time and sampling the average absorption line 
profile with the FWHM of » 15 km s '. Using simulations, we 
found that for small inclinations of the planetary orbit, where the 
planets appear only slightly illuminated, the method is unable to 
detect Jupiter-size objects with very high albedos. 



4.5. Step 4: Determination of the confidence level 

Once the best model M[Kp,e(A)] has been evaluated, we de- 
termine the confidence level of the x^ minimum by applying 
the bootstrap randomisation method (e.g. Kiirster et al. 1997). 
Retaining the orbital phases, we randomly redistribute the ob- 
served spectra amongst the phases, thereby creating A^ different 
data sets. Any signal present in the original data is now scram- 
bled in these artificial data sets. For all these randomised data 
sets, we again evaluate the model for the two free parameters, 
and locate the best fit with its specific x^ minimum. We set m 
to be the number of best-fit models to the A^ randomised data 
sets that have a minimum ;^'^ less or equal than the minimum ;^'^ 
found for the original data set. The confidence level can then be 
estimated by x; 1 - m/N. 



5. Results and discussion 

Applying the data synthesis method to the HD 75289A data, we 
adopted the following approximations to the atmospheric mod- 
els by Sudarsky et al. (2000). 

(i) We adopted a grey-albedo model to resemble the Class 
V model, which describes the atmospheres of hot Jupiters with 
temperatures > 1500 K. As can be seen in Figure [1] this was a 
valid approximation for our observed wavelength range (402 to 
522 nm). 

(ii) We also considered the two Class IV models for hot 
Jupiters with temperatures around 1300 K. The irradiated model 
predicts a very low reflectivity due to the lack of silicates in 
the upper layers of the planetary atmosphere, while the iso- 
lated model predicts high reflectivity due to the existence of sil- 
icate clouds in the upper decks of the planetary atmosphere. We 
found that both Class IV models show a similar trend for the 
albedo function. Therefore we considered only that trend in fur- 
ther modeling. 



5. 1 . Grey albedo model (corresponding to Class V) 

Our data analysis, using a grey albedo, revealed a;^f^-minimum 
at a planet-to-star flux ratio e(A) = 2.9 x 10^^, and an RV semi- 
amplitude Kp - 121 km s ' corresponding to an orbital inclina- 
tion of / - 54°. Figure |2] (upper panel) shows ax^ contour map. 
However, using bootstrap randomisation with 3000 trial data sets 
we found that this ;^'^ -minimum was uncertain with a false-alarm 
probability (FAP) of 19 %; we do not therefore consider this as 
a detection of reflected fight from the planet. 



5.2. Class IV model 

The analysis using the Class IV albedo model did not yield any 
evidence for reflected light from the hot Jupiter HD 75289Ab. 
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Fig. 2. Contour maps of x^ normalised to x^ for the model 
parameters Kp and e{A) showing the result of the analysis of 
the HD 75289A data for two different atmospheric models. For 
better visualisation the x^ contour levels follow the sequence 
xl=x'-- +0.1;f-2 e*-i forA:>0. 

^ K ^ mill ^ min 

Upper panel: For the planetary model, we assumed a Venus-like 
phase function and a grey albedo. The minimum ;^'^ was found at 
/Tp = 121 km s'' and e{X) = 2.9 x 10 ^ An analysis with 3000 
trials revealed that this feature is uncertain with a false alarm 
probability of 19%. 

Lower panel: For the planetary model, we assumed a Venus-like 
phase function and a class IV albedo function. The minimum 
X^ was found at /Tp = 119 kms^' and e(^) = 3.3 x 10"^ An 
analysis with 3000 trials revealed that this feature is uncertain 
with a false alarm probability of 24%. We note that the corre- 
sponding geometric albedo p is shown under the assumption that 
^p - 1.2/?jup (the right-hand y-axis). 



We found that the ;t'^-minimum at a planet-to-star flux ratio 
e{A) = 3.3 X 10"^ and an RV semi-ampUtude /Tp = 119 km s"', 
which corresponds to an orbital inclination / - 53°. Figure |2] 
(lower panel) shows a contour map of that feature. From boot- 
strap randomisation with 3000 trials we found that this feature 
was uncertain with a FAP of 24%. 



5.3. Upper limits 

Obviously, the amount of data for when the planet was observed 
during its bright phases was insufficient to measure the plan- 
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Fig. 3. Contour map showing confidence levels for the upper lim- 
its to the planet-to-star flux ratio e{A) as a function of to the RV 
semi-amplitude Kp of the planet (lower x-axis), or orbital in- 
clination / (upper X-axis). From top to bottom, these levels are: 
99.9, 99.0 and 90.0 % confidence. For the planetary model, we 
assumed a Venus-like phase function as well as a Class IV albedo 
function (solid lines) and a grey albedo (dashed lines). The cor- 
responding geometric albedo p is shown under the assumption 
that Rp - 1.2 R]ap (the right-hand y-axis). The region above the 
thick horizontal line is excluded because the albedo of the plan- 
etary atmosphere would exceed 1 . 



etary signal. The question is now: at what planet-to-star flux 
ratio e{A) would we be able to detect significantly the plane- 
tary signal using the present data? To this end, we constructed 
data sets by adding an artificial planetary signal to the original 
data. To produce the planetary signal, we explored the range of 
Kp e [40, 180]kms-i ande(A) 6 [2x10^^ 2.5x10-"^]. For each 
data set, the simulated planetary signal was a copy of the object 
spectrum, but shifted by Vp(Kp, (p) and scaled down by the factor 
/(0, /, A), which includes the albedo model of the planetary at- 
mosphere (either grey albedo or Class IV model). For each pair 
of the parameters .^p and e{A), we determined the confidence 
level of the best-fit model to the artificial data via bootstrap ran- 
domisation. 

This enables place upper limits to be established to the 
planet-to-star flux ratio e{A) for different confidence levels, as 
a function of the RV semi-amplitude Kp. Figure [3] shows that 
the upper limits to the planet-to-star flux ratio decrease with in- 
creasing orbital inclination, which is a direct consequence of the 
illumination geometry. As can be seen from Figure [3] (dashed 
lines), we have the highest sensitivity for detection of the plane- 
tary signal at high orbital inclinations with the formal numerical 
optimum found at an RV semi-amplitude Kp - 129 km s"' cor- 
responding to / = 60°. At this inclination, the 99.9 % confidence 
upper limit to the planet-to-star flux ratio for the grey albedo 
model is e = 6.7 x 10"^ , while it is e = 8.3 x 10"^ for the Class 
IV model. In Figure |3] we clearly see that the upper limits es- 
tablished by adopting the grey albedo model are deeper than the 
ones found with the Class IV model. A plausible explanation for 
this seems to be that in the grey albedo model, all the absorp- 
tion lines have the same weight. On the other hand, in the Class 
IV model, the absorption lines in the red part of the observed 
spectrum have a lower weight than those in the blue part. Hence, 
fewer lines are effectively used for the analysis. 
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Table 3. Estimates for the observing times required for UVES 
to collect high-S/N data sufficient to detect, with 99.9% con- 
fidence, the reflected light from HD 75289Ab with a specific 
planet-to-star flux ratio e. Here we assume a Class V model for 
high-reflective case and the isolated Class IV albedo shape for 
the fainter case (cf. Fig. 1). The corresponding values of the 
albedo values p are based on the assumption that the planetary 



radius is Rp - 


1.2/?jup. 






e 


p N (spectra) 


time (h) 


albedo model 


5.7 X 10' 
2.9 X 10-^ 


0.4 324 
0.2 1460 


19.5 
87.6 


Class V 
isolated Class IV 



Assuming a planetary radius Rj, - 1.2 Rj^p, we find that the 
99.9 % upper limit to the geometric albedo is p < 0.46 for the 
grey albedo model, and p < 0.57 for the Class IV model, at 
a wavelength of A = 445 nm, which corresponds to the center 
of gravity of the albedo function. For comparison, at this wave- 
length Jupiter's albedo is p - 0.43 (Karkoschka 1994). We note 
that the planetary radius of /?p = 1 .2 Rj^p is the average value of 
the five transiting hot Jupiters listed in Section 2.3. 

In the case that the planetary radius of HD 75289Ab is simi- 
lar to that of the transiting hot Jupiter HAT-P- 1 b (/?p = 1 .36 Rjup, 
Bakos et al. 2007), we find the geoemetric albedo to be /? < 0.36 
and p < 0.45, respectively for the grey albedo model and the 
Class IV model. This would enable us to rule out the Class V 
albedo model. 



6. Conclusions and outlook 

We reanalysed 684 high-precision spectra of the star 
HD 75289A, obtained during four nights by another re- 
search group in 2003, in an attempt to detect starlight reflected 
from the known hot Jupiter. 

(i) We found that the upper limits to the planet-to-star ra- 
tios established by Leigh et al. (2003a) were based on erroneous 
orbital phase information, and therefore needed to be corrected. 

(ii) Our reanalysis using the correct orbital phases produced 
a non-detection of the reflected light of the planet. We therefore 
placed upper limits on the planet-to-star flux ratios for different 
possible orbital inclinations and confidence levels. These upper 
limits, however, were s; 60 % higher than the ones placed by 
Leigh since the planet was not observed at its best visibility. 

(iii) We furthermore determined values of the geometric 
albedo that was 2.9 times higher than the values established by 
Leigh et al. (2003a), since, due to limited information on plane- 
tary transits available in 2003, these authors assumed a planetary 
radius that was larger by 33 % than our value. 

(iv) The previous upper limit to the geometrical albedo of 
0.12 established by Leigh et al. (2003) was by far the deepest 
one determined for hot Jupiters up to now. These authors sug- 
gested that a Class V atmospheric model can be ruled out for 
HD 75289Ab. For the value of the planetary radius adopted by 
us, Rp - 1.2 Rjup, we find that the amount of data was insuf- 
ficient to rule out specific atmospheric models. However, if the 
radius of HD 75289Ab is similar to that of the transiting planet 
HAT-P-1 b (Rp = 1.36 R},,p, Bakos et al. 2007), we can rule out 
the Class V albedo model. 

To detect reflected starlight from the hot Jupiter 
HD 75289Ab, more high-precision data would be needed. 
Table [3] lists estimates of observing times required to detect, 
with 99.9% confidence, the reflected light of a planet with a 



Class V model and an isolated Class IV model, respectively. 
These calculations are based on the following assumptions: the 
fictitious observations are scheduled optimally with respect to 
phase (0 - 0.30 to 0.45 and 0.55 to 0.70). These observations 
are carried out by using the blue arm of UVES at the VLT/UT2 
in Chile. The wavelength range is /I = 400.3 to 527.0 nm, and 
the exposure time is chosen in such a way that a maximum 
count rate of 40 000 ADU/px is achieved (cf. Section 3). 
Table [3] demonstrates that we would require about 20 hours of 
observations to measure the reflected light of HD 75289Ab in 
the case that it had a Class V albedo function and a planetary 
radius of /?p = 1 .2 Rji,p. However, if HD 75289Ab had an albedo 
function similar to that of the isolated Class IV model, it would 
already take more than 87 hours of observing time to detect the 
reflected light. 

Because of the large amount of data required, we conclude 
that the star HD 75289A is too faint for the search for re- 
flected starlight. Instead, the search for thermal emission of 
HD 75289Ab in the near-infrared appears to be promising, since 
the planet-to-star flux ratio increases dramatically in the infrared 
(Barnes et al. 2007). 
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